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ABSTRACT 

Shown is that contrary to common intuition, even an arbitrarily weak at- 
tenuating mechanism is sufficient to make the background sky quite dark inde- 
pendently of the size of the universe and the Hubble expansion. Further shown 
is that such an attenuation already exists in the wave nature of light due to 
entrapment and diffusion from successive diffractions. This is a fundamentally 
new mechanism to physics, as illustrated by application to the solar neutrino 
attenuation, galactic dark matter and gamma ray bursts problems. It not only 
provides a big bang-like cutoff, but also appears to explain the appearance of 
primeval, metal-deficient galaxies at high redshifts, without deviating from the 
Olbers' premise of an infinite universe. 



1. Introduction 

This historical paradox highlighted by Olbers almost two centuries ago (Olbers 1826) 
holds that the night sky should have been as bright as the sun's disk because we should en- 
counter rays from a stellar surface in any direction we look. More particularly, the argument 
holds that the stellar light must cover an increasing area oc r 2 as its distance r from its source 
increases, its brightness J diminishing as Jo/r 2 , where Jo is the stellar surface luminosity. 
This is the well known inverse-square law for light, and takes only the geometrical spread 
in 3-dimensional space into account. According to the paradox, as we look at regions of the 
sky where the stars are very far, they should each appear dimmer, but a given solid angle 
would cover oc r 2 stars, assuming a uniformly populated infinite universe. The luminosity of 
the sky should therefore be cx (J /r 2 ) x r 2 = J , meaning that the background sky should 
be so bright that the stars should be indistinguishable against it. I present two fundamental 
results, first, that an infinitesimally small attenuation a > 0, barely enough to change the 
propagation law to the form J e~ ar /r 2 , suffices to solve the paradox, and second, that such 
an attenuation happens to be inherent in the wave nature of light. 
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The results are unintuitive because almost any kind of attenuation leads to this form, 
and absorption and scattering by dust have been considered inadequate in the past. Harrison 
argued (Raychaudhuri 1979) that the dust would eventually attain thermal equilibrium with 
the stars and effectively stop absorbing more energy. Since the radiation too would even- 
tually reach equilibrium, scattering of itself cannot solve the paradox either. The standard 
model offers three plausible solutions that the Hubble flow causes the light to lose energy 
well in excess of the inverse-square attenuation, that the universe is too young for thermal 
equilibrium, and that the universe is as such finite. Wesson has shown that the first would 
actually contribute less (Wesson et al. 1987), so finiteness of the universe, in both age and 
extent, is currently believed to be the reason for the darkness of the night sky. I shall show 
that a similar cutoff occurs because of the inherent attenuation, that it necessarily leads 
to a further mechanism of spectral modification that would make the most distant galax- 
ies appear primeval, which is again known and currently attributed to the big bang. In 
both cases, I exploit the smallness of the mechanisms needed for the respective effects. The 
approach is only possible because of the exponential factor in the attenuated propagation 
law, but more importantly, it demonstrates that the small orders routinely left out in the 
approximations applicable on earth could be significant physics on the cosmological scale. 
To further emphasise the error in letting approximations dictate our reasoning, I shall show 
that an even smaller order of attenuation due to the same mechanism would explain the 
missing solar neutrino flux, currently attributed to neutrino oscillation, and could mean new 
insight of a fundamental kind in the physics of weak interactions. 

Unlike the dust effects, the attenuation of present concern depends only on the presence 
and not the thermal state or other properties of matter, and is therefore immune to thermal 
equilibrium. It results from a diffuse entrapment of radiation due to successive diffraction 
and gravitational deflections that continually turn a portion of the wavefront. While the 
occurrence of successive deflections is known, for example, in Laue diffraction theory, its 
implications to astrophysics have not been examined in previous treatments, such as in the 
context of extinction by dust (Spitzer Jr 1978, pl49-153), possibly because one ordinarily 
thinks of diffraction as carrying wave energy around obstructions, increasing rather than 
diminishing the net power flow. This would be the case in a hypothetical universe where 
the sources are assumed to be behind a plane of diffracting obstructions, but the notion 
does not really extend to the three-dimensional universe involved in the paradox, where the 
stars themselves are the obvious obstructions to each other's light. The enhancing property 
turns out to be intransitive because the successive deflections can then keep a fraction of the 
radiation from ever reaching its original destination. 

Another reason why this result was unobvious is that Fraunhofer's approximation is 
invariably assumed because of the immense distances involved. A solution to Olbers' paradox 
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then appears to be ruled out because the loss of the direct rays from a star due to en route 
diffraction by an angle say 9, would be made up by the rays from another star behind the 
diffractor at an angle —6. While the argument is somewhat weaker than enhancement, it 
reveals the error in our past intuition, because the probability p(r s ) of finding a compensating 
star at the same or less distance r s behind the diffractor depends on r 2 /r 2 , which is certainly 
less than unity. The approximation assumes that both the sources and the observer are 
infinitely far from the diffracting object, i.e. r s — > oo and r — > oo, in which case the ratio 
does not matter. These conditions are, however, applicable only in the vicinity of a given 
diffractor, and cannot be legitimately applied when the obstructing objects are distributed 
over the same scale of distances as the sources. In the presence of an attenuation a, the 
distance r s matters because the compensating source becomes more likely to be farther by 
the triangle theorem, and therefore likely to be dimmer. 

In both the real universe and Olbers' scenario, therefore, the light from a distance 
source located on a geometrically unobstructed straight line from us does get diminished by 
diffraction due to obstructions lying off the straight line. This loss would be compensated, as 
in the Fraunhofer case, if enough diffracted light from elsewhere could rejoin the straight line 
path. Traditional wisdom suggests that the repeated deflections be treated as as a random 
walk leading to a slow diffusion of the photons, which would not yield a net reduction 
of the average luminosity. However, there are problems with this view, because while it 
is reasonable to assume that the diffracting objects are randomly distributed, their gross 
motions are not random and are quite slow in relation to the interstellar distances. 

What we have is an essentially static pattern lacking the temporal randomisation of 
direction needed to qualify it as random walk. More particularly, all efforts to simulate 
thermalisation with fixed dynamical models have consistently led to persistent oscillatory 
states, since the very first attempt in 1953 by Fermi, Pasta and Ulam (Fermi et al. 1965; 
Fillipov et al. 1998), showing that mere complexity of dynamical structure is not enough 
for assuming diffusion. Moreover, any static pattern necessarily contains circulations, which 
might not only explain the FPU problem, but in our case, trap some of the light virtually 
forever. Harrison's argument cannot be applied to such states because the circulations would 
be centrally dependent on the individual sources, and each of which presumably has a finite 
lifetime even in Olbers' scenario. We do expect most of the "trapped" energy to eventually 
diffuse out, but we have no basis to assume that all of it will. Rather, we can expect a portion 
of the energy to get absorbed or turn into matter, given that the attenuated propagation 
law is already characteristic of the Klein-Gordon equation (V 2 + a 2 d 2 / 'dt 2 ) ip — 0. The latter 
is simply the quantum version of the relativistic argument that radiation when retarded to 
effective speeds less than c should exhibit rest mass, and relates to our treatment of the solar 
neutrino problem. We can thus be certain of a net attenuation a > 0, and this suffices, as 
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shown next, to solve the paradox. 



2. Solution of the paradox 



As stated, our key argument is that any attenuation whatsoever, so long as it operates 
on the large scale, solves the paradox independently of the standard model. We seek an 
attenuation a (dB m _1 ), such that the propagation law for light changes from r~ 2 to r~ 2 e~ ar , 
which reduces the background brightness of the sky to 



While the integral is well known, the solution is not immediately obvious, largely because 
the value superficially resembles Jo, making it look as if we need a very large a, of the order 
of at least 130 dB = 10 13 (Roach and Gordon 1973, p. 24-25) to get a dark background sky, 
and hitherto seemed impossible without the big bang theory. 

The resemblence is misleading because the unattenuated J differs in dimensions from 
Jo/a, as the latter has the dimensions of luminosity x distance. For legitimate comparison, 
we must express J in exactly same dimensions, hence in the statement of the paradox, the 
"Olbers luminosity" J, which one may informally think of as the brightness of the sun's disk, 
corresponds not to a — 1 but to a — 0, i.e. 



which is infinity. Conversely, when we calibrate with respect to the sun's disk, the attenuated 
background sky should be infinitesimally dim, for any o > 0. To appreciate why this should 
be so, consider how bright the background needs to be in order to match the sun's disk. 
In Olbers' argument, whichever direction we look in, our line of sight must meet a stellar 
surface and at any finite angular resolution 9, the brightness should correspond to the number 
of stars included within the solid angle 9. Accordingly, the observed brightness would be 
J rj = Jo eT ar along that direction, and that set of stars then needs to be e ar times brighter 
than the sun in order to match Jo. Clearly, it does not matter if a is very small; as long as 
it is nonzero, we have an effective cutoff of the observable universe at 




(1) 




(2) 



r n = a 1 log( J /J n ) 



(3) 



for a finite a, where J n represents the background noise in the measuring process. By eq. (3), 
the standard model cutoff, corresponding of an age of the universe of ~ 15 Gy, is equivalent 
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to 



a ~ 130 dB/15 Gy = 9 x 10" 25 dB m -1 . 
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This is far less than what one might naively expect from the form of J a (eq. 1), and amounts 
to a mere 8 x 1CT 9 dB per light-year. As remarked in the introduction, it is the smallness of 
the attenuation needed to explain the apparent big bang cutoff that makes it impossible to 
rule it out on the basis of terrestrial physics. Furthermore, eq. (4) is as yet an upper bound, 
because we ignored the diffraction from nonluminous bodies as well as the gravitation of 
these and the visible objects, which contribute substantially to the deflections. We also 
ignored the dust extinction in our own neighbourhood (Roach and Gordon 1973, ch.4), to 
which Harrison's argument again does not apply, and the impact of the Hubble flow, which 
together account for a good part of the 130 dB. Note that Harrison's argument remains valid 
for dust on the large scale, and the subtlety that overcomes it for the diffractive scattering, 
described in the next section, is the coherence and dependence of the diffuse circulatory 
states on their respective sources. 



Eq. (2) establishes our principal point that only an extremely small attenuation is needed 
to reproduce the big bang cutoff. It remains to be shown that such an attenuation is indeed 
possible and likely from successive diffraction. As explained in §1, we are concerned with 
multiple interstellar hops that are each much larger than stellar diameters, so that the 
Fraunhofer theory can be applied to the individual diffractions. We may more particularly 
treat the stars as point objects, representing the diffraction by an angular spreading function 
f(4>, 5), f > V {(f), 9}, applicable to a parallel beam of incident light, giving 



for the diffracted light, being the azimuthal angle around the beam axis and S, the angular 
spread from the axis. We shall now examine the special treatment necessary to account for 
successive diffractions. 

The spread function / is rather like the differential scattering function a(Q) of the 
Rutherford model (Goldstein 1980, §3-20), but several important differences must be noted. 
Firstly, we are concerned with continuous wavef unctions, not particles, so that the "scatter- 
ing" itself is not at all probabilistic; the only probability inherent in the model is the stellar 
distribution. Secondly, the scattering function / depends on the X/D ratio, where A is of 



3. Scattering approximation 




(5) 
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course the wavelength and D, the mean stellar diameter. This is treated in more detail in 
terms of Fresnel-Kirchhoff theory in Appendix A. More importantly, the result of interest is 
the net attenuation and not a cross-section of interaction. It is common to use a for both, 
but the second notion is not of concern here. We shall now formalise the cumulative forward 
loss and the nonzero "backscatter" occurring at each encounter. 

The second of eqs. (5) represents the conservation of energy at each diffractive encounter, 
since it means that 

/ J(<j>,6)d<j>dS = Jo. (6) 

r =0 J 5=0 

We start by applying / to a bundle of rays incident on a first star. Since the rays spread as 
if from a point source, they acquire a l/47rr 2 spreading loss, even if we had started with a 
parallel bundle, which would be equivalent to a planar wavefront. A second encounter after 
a distance r\ therefore yields 



J{M) = iTi f M9-8)f((l>,8)d6 
4vrr x J s=0 



(7) 



(8) 



(9) 



where g\ denotes the first cumulative integral 

gi(<P,9) = ±- I f(M-8 x )f(<l>M)dh. 
4tt J 5l=0 

We thus obtain the recursive set of integrals 

g = f and 

9^ d ) = ^J f{<t>,0-S)9n-i{4>,5)d5 for n = 1,2,... 
describing n successive encounters as 

n 

J(0, 6) = Jo g n (<f>, 6) J] - » J g n (<t>, S) f~ 2n (10) 

j=l r J 

where r 3 - are intervening distances, and f, the mean distance between stars. Integrating eq. 
(10) over 0, and summing over the contributions from one or more encounters, we obtain 
the total "gain" at 9 from the initial beam as 

n(a\ #o( g ) i n i9i{9) n 2 g 2 (9) 

W = -=2~ + ~A + ^ + - 

L (11) 



where g n (9) = g n ( 
Jo 
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Here, rij > 1 denote the mean number of parallel paths corresponding to the range of 
for j encounters, and gj(9) are the per-path contributions. It should be at least intuitively 
clear that rij would be an increasing combinatorial function of j, offsetting the increasing 
geometrical attenuation from the denominator (f -2 -'). 

Eq. (11) reveals an interesting property of the successive diffractions: that the Qj(9) 
increase with j for precisely the reason that non-paraxial angles are ordinarily ignored in 
terrestrial optics - at large j's, small incremental angles 5 become significant in the integrand, 
so that 



as the integral converges to the zero-th order beam. Eq. (12) represents the observation that 
a succession of small diffraction angles adds up to a large total deflection with significant 
amplitude, since the component deflections are paraxial. Three encounters each of 7r/3, 
for example, suffice to make a back contribution, and the direct backscatter g (9 ps n) 
suffers no 1/f 2 attenuation as it involves no intermediate hops at all. As a result, nonzero 
total reflection is generally guaranteed and the increasing rij also partly compensate for 
the 1/f 2 - 7 loss. 5'o( 7r ) is generally ignored in terrestrial optics because the cos 2 (6?) factor is 
more pronounced over small distances; this too is inappropriate over the cosmological scale 
of distances. As stated in §1, motions of the surface of the diffracting star, although many 
orders larger than the wavelength of light, is not significant in the present context, principally 
because / represents angular distribution of the radiant power and is unaffected by phase 
fluctuations. 

The presence of a net backscatter g(n) means that light from sources in the observer's 
own neighbourhood will tend to lighten the sky. This is an expected effect of scattering due 
to dust, for example: Roach and Gordon estimate that for an observer near the centre of 
our galaxy, the night sky should be relatively opaque (Roach and Gordon 1973). This does 
not contradict either the paradox or Harrison's argument, but the diffractive backscatter 
potentially does because it would be spread over cosmological distances. Our result remains 
intact because the backscattered light would itself be again subject to repeated diffractions 
and suffer the same attenuation. The totality of "local" sources we need to consider would 
cumulatively increase as r 3 , but gets overcome by the exponential attenuation e~ ur within a 
limited r. 




(12) 
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4. Applications 

As mentioned in §1, our mechanism has several interesting applications, beginning with 
the manifestation of rest mass of the photons because of the Klein-Gordon form of the 
attenuated propagation law. This suggests that our trapped circulatory photon states could 
be responsible for a portion of the dark matter especially if the universe were very old. A 
related observation is that the density of circulatory states would be proportional to the 
density of particulate matter, because the deflections would become available at smaller 
distances in dense neighbourhoods. We would expect the attenuation to be greater within 
galaxies than in the sparser regions between them, and the resulting "optical dark matter" 
distribution to be consistent with the overall galactic dark matter indicated by their rotation 
profiles, as the circulatory states would presumably add up with radial distance from the 
galactic centre. At present, we do not know how to estimate this optical component of dark 
matter. 

The attenuation might also explain the intensity of gamma ray bursts from galactic re- 
gions. For a given aperture between the stars and other objects surrounding the line of sight, 
the basic principles of diffraction dictate that gamma rays would suffer far less diffractive 
loss than visible light because of their shorter wavelengths. The observed gamma ray burst 
intensities should therefore be much more than we might expect from the visible brightness. 
While the reasoning appears to be in the right direction, it remains to be substantiated. 

Another high energy scenario consistent with our theory is the known attenuation of the 
solar neutrino flux, which, as mentioned in §1, is currently attributed to quantum oscillation. 
Given that the mass of the sun is M Q 2 x 10 30 kg and the average mass of a nucleon, 
1.67 x 1CT 27 kg, we estimate there are 1.2 x 10 57 nucleons within the sun's volume ~ 1.4 x 
10 27 m 3 , which allows roughly 1.2 x 10~ 30 m 3 per nucleon. This means that a neutrino would 
encounter a nucleon every 1.06 x 10 -10 m, or up to 6.6 x 10 18 times from the centre to the 
surface. This mean free path between nucleons is the deBroglie wavelength for a particle of 
about 12 keV, so some diffraction appears to be inevitable. The observed attenuation, about 
two-thirds, must result from O(10 18 ) compoundings, assuming, as an approximation, that all 
neutrinos originate from near the centre; this yields a exp[lCT 18 log(2/3)] exp(4 x 1CT 19 ) 
or 10~ 19 dB per encounter, as the diffractive attenuation needed to explain the observed loss. 
Once again, it is the smallness of the necessary attenuation that makes it plausible as the 
cause, and we speculate that this could be the wave-theoretic reason for the breaking of 
symmetry that leads to weak interactions. 

Our last application is to the appearance of primeval galaxies at high redshift factors z, 
and comes from observing that the selective absorption of certain frequencies by matter close 
to the line of sight should trigger diffractive loss a' ps a at the absorbed frequencies. This 
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should result in a subtractive spectral modification of the light from more distant matter, 
and the spectral subtraction by our own galactic neighbourhood, which is well evolved and 
metal-rich, should make more distant galaxies appear metal-deficient and primeval. Only 
an extremely small a' is needed, once again, to explain the almost total metal-deficiency of 
the primeval galaxies. As a 1 must result from selective absorption by obstructions that must 
be transparent at other frequencies, i.e. from a different set of obstructions, a' cannot be 
identical to a. But it would be of the same order, leading to a similar cutoff r' n r n (eq. 3), 
because denser regions of dust that invariably surround opaque bodies are likely causes of 
selective absorption. Importantly, we may argue that the most distant galaxies can appear 
mature only if there were a uniform mix of young and old galaxies all the way right up to our 
own galaxy, but this would contradict the standard model. This does not mean unexpected 
support for the big bang theory, however, because we could merely be in a relatively small 
region of the Olbers' universe, of radius at least r n , that formed at or before the big bang 
indicated by the Hubble redshift. 

I thank Shai Ronen, Gyan Bhanot, A Joseph Hoane and most of all, Bruce Elmegreen, 
for valuable discussions in the context. 

A. Point diffraction 

The usual treatment for diffraction due to a small region concerns apertures in an 
opaque screen, but our stars are in effect obstructions in an otherwise empty space. If 
U s represents the diffracted amplitude due to a star, and U a , that due to an aperture of 
the same size, we have by Babinet's principle that U s = Uq — U a , where Uq represents the 
unobstructed incident wavefront. In holography, for instance, we can use this theorem to 
compute the resulting interference between the direct rays (Uq) and the diffracted rays (U s ) 
from a point obstruction. In the present context, however, we are not interested in the 
local interference patterns, but only in the direction of the power flow, so the interference 
between direct and diffracted wavefronts is irrelevant. We do sum over rays going in the 
same direction 8 in eq. (11), but f ^> A, and in any case, there is enormous thermal jitter 
from our stellar obstructions, which would randomise the phase between the summed paths, 
hence interference effects can again be ignored, even if our initial rays happened to be quite 
coherent. 

We do need to consider phase within narrow bundles of rays, however, in order to obtain 
the wavelength dependence of the diffraction. Under the far-field assumption f ^> D, D being 
the mean diameter of our stellar obstructions, the incident wavefronts before diffraction 
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would be not only almost planar, but likely to be temporally coherent across a sizeable 
fraction of the stellar surface. We shall limit ourselves to pure Fourier components, avoiding 
questions of incoherence, because issues of decoherence and wavepacket dispersal can be 
considered in terms of these components. With these assumptions, we may apply Kirchhoff 's 
condition that the conditions in a stellar diffracting region A will not be much affected by 
the surrounding matter, which yields the usual Fresnel-Kirchhoff diffraction formula (Born 
and Wolf 1959, §8.3.2) 

j f p ik(r+s) 

U ^ = " *T / t cos ( f ) " cos (*)] dS ' ( A1 ) 

2A J A rs 

where s and r are distances to the source and to the point of observation, respectively. We 

cannot, however, afford to replace the cosines with a single cos(#) factor as in the traditional 

theory, since this would be proper only for paraxial rays. To correctly compute U at large 

diffraction angles 5, we need to retain the two-cosine form in our Fraunhofer approximation, 

obtaining 



= _ i [008(f) - COs(g)] j (r+-) ^ 

_ _ iu [cos(f) - cos(g)] e ik{r+s ^ r j k{ z +r)) dS 
27rc rs J A 



(A2) 



£ and i] being the direction cosines as in prior theory. A remaining problem is, of course, 
that the above formulae are meant for diffraction from an aperture, and not an opaque 
disk presented by each star, but Babinet's principle assures us that the resulting / will be 
identical. Assuming our stars to be generally circular, we find that the integrand would have 
the form 2Ji(ka9)/ka9, where k = 2n/\ and a is the radius of the aperture, yielding 

"2Ji(Jfea0)" " 
(2ttc) 2 [ ka6 

since the (rs) 2 factor is separately accounted for by f in §3. 



/(fl)~|^)f| 2 = ^ [008( ^- 008(5)11 



(A3) 



As a result, f(9) is non-zero almost everywhere, over the entire range of 6, which in- 
cludes 7r. This "backscatter" is a reflection from space, due to the local perturbation in the 
impedance of space caused by the off-axis diffracting star, and is definitely nonzero for any 
nonzero solid angle around 9 = ir. The angular spread is more strongly governed by the 
phase factor than by u>, so that a given gap between nearer obstructions will be more open 
to gamma rays than to visible light, as remarked in §1. 
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